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Abstract:  We  apply  the  equivalent  theory  to  orthorhombic  anisotropic  materials  and  provide 
a  general  unit-cell  design  criterion  for  achieving  a  length-independent  retrieval  of  the  effective 
material  parameters  from  a  single  layer  of  unit  cells.  We  introduce  a  graphical  retrieval  method 
and  phase  unwrapping  techniques.  The  graphical  method  utilizes  the  linear  regression  technique. 
Our  method  can  reduce  the  uncertainty  of  experimental  measurements  and  the  ambiguity  of 
phase  unwrapping.  Moreover,  the  graphical  method  can  simultaneously  determine  the  bulk 
values  of  the  six  effective  material  parameters,  permittivity  and  permeability  tensors,  from  a 
single  layer  of  unit  cells. 
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1.  Introduction 


Metamaterials  (MMs)  are  artificial  materials  engineered  to  achieve  unusual  electromagnetic  (EM)  properties  that  are 
not  normally  found  in  nature ifTl-l^.  A  comprehensive  review  can  be  found  in  the  booklU.  Distinguished  from  photonic 
crystals,  metamaterials  are  made  from  periodical  structures  with  unit  cells  much  smaller  than  the  wavelength  of  light. 
A  general  method  to  construct  large  area  3D  MMs  is  layer-by-layer  fabrication  technique  111- HI.  In  stacked  MMs, 
interaction  between  adjacent  layers  makes  it  difficult  to  extract  bulk  material  parameters  from  a  single  unit-cell  layer. 
It  has  been  found  that  the  retrieved  effective  metamaterial  parameters  are  often  dependent  on  the  number  of  unit  cells 
along  the  propagation  direction ||5l|9l.  It  appears  that  there  is  no  clear  methodology  on  how  to  accurately  predict  the  bulk 
values  of  the  effective  permittivities  and  permeabilies  through  a  single  layer  of  unit  cells.  Currently  there  is  no  simple  and 
effective  way  to  resolve  phase  ambiguity  in  determining  the  phase  of  the  transmission  and  reflection  of  electromagnetic 
fields  because  of  phase  wrapping.  Although  the  phase  ambiguity  is  a  common  issue  for  the  parameter  retrieval  of  the 
general  composite  materials,  this  issue  becomes  more  significant  for  metamaterials  where  typically  resonances,  positive 
refractive  index,  and  negative  refractive  index  are  all  present  in  the  same  frequency  band.  In  this  paper,  we  apply 
Herpin’s  equivalent  theoremlflOl  to  orthorhombic  anisotropic  media  and  provide  a  simple  way  to  accurately  predict 
the  effective  bulk  material  parameters  from  a  single  layer  of  unit  cells.  We  introduce  a  graphical  retrieval  method  and 
phase  unwrapping  techniques,  which  can  simultaneously  determine  the  six  material  parameters,  the  permittivity  and 
permeability  tensors,  from  one  unit  cell. 

2.  Equivalent  theory 

Layer-by-layer  fabrication  method  renders  metamaterials  intrinsically  anisotropic.  In  MMs  one  unit-cell  layer  is  usually 
composed  of  several  sub-layers  of  different  materials  or  nanostructures.  We  limit  our  discussion  in  nonchiral  materials. 
To  illustrate  the  key  point  we  adopt  multilayer  method  that  has  been  used  to  model  MMs  by  many  groups  ifTTI-lfTTl.  We 
further  assume  that  the  principal  axes  of  the  sub-layers  are  parallel  in  each  direction.  In  the  principal  coordinate  system, 
the  permittivity  and  permeability  tensors  are  given  by 


Enx 

0 

/  l^nx 

0 

e„  =  0 

^ny 

0  , 

A«  =  0 

Mny 

0  , 

(1) 

VO 

0 

VO 

0 

P-nzJ 

where  n=  1,2,  •  •  • .  The  scalar  terms  Enj  and  {j  =x,y,z)  are  complex.  Consider  a  monochromatic  wave  of  frequency 
0)  with  time  dependence  exp{—i(0t)  propagates  inside  the  orthorhombic  anisotropic  materials.  In  each  layer  we  have 

Wx{^-^-WxE)  =  kl{e„-E), 

V  X  -  V  xH)  =  kl{fln-H) , 

where  kQ  =  (O /c.  If  the  plane  of  incidence  is  one  of  the  crystal  planes,  the  TE  and  TM  polarizations  are  decoupled. 


Fig.  1.  A  schematic  shows  how  to  construct  a  symmetric  unit  cell,  which  is  composed  of 
two  original  asymmetric  unit  cells. 

According  to  Herpin’s  equivalent  theoremlflOl.  every  general  multilayer  is  equivalent  to  a  two-homogeneous-layer 
system  and  every  symmetric  multilayer  is  equivalent  to  a  single  homogeneous  layer,  characterized  by  an  equivalent 


index  and  equivalent  thicknessifTSll.  In  metamaterials,  a  single  unit-cell  layer  can  be  considered  as  a  multilayer  system. 
So  it  is  equivalent  to  a  two-homogeneous-layer  system,  denoted  as  AB.  We  can  then  construct  a  symmetric  unit  cell  by 
cascading  two  unit  cells  as  ABBA,  which  is  equivalent  to  a  single  homogeneous  layer  according  to  Herpin’s  theorem. 
This  process  is  illustrated  in  Fig.  [T]  Applying  this  methodology,  a  symmetric  unit-cell  layer  can  often  be  constructed 
regardless  the  number  of  sub-layers  and  complexity  of  each  sub-layer  in  the  original  unit  cell.  The  permittivities  and 
permeabilities  retrieved  from  the  symmetric  unit  cell,  which  is  composed  of  two  original  asymmetric  unit  cells,  will 
represent  the  bulk  material  parameters.  Thus,  a  length-independent  description  can  be  achieved. 


To  make  it  more  clear,  let  M  represents  the  characteristic  matrix  of  one  symmetric  unit  cell.  According  to  Herpin’s 
equivalent  theorem,  it  can  be  replaced  by  an  equivalent  single  layer.  Assume  the  x-z  plane  is  the  plane  of  incidence.  For 
TM  mode,  i.e.  H  =  (p,Hy,0)  and  E  =  {Ex,0,Ey),  we  have 


Mil  =M22  =  cosV4, 

sin  We  (3) 

=  M2i=-iXsm\t/e, 

where  1/4  is  the  equivalent  phase  thickness  of  the  symmetric  unit  cell;  is  the  equivalent  impedance.  If  the  material 
contains  A-layer  symmetric  unit  cells,  the  characteristic  matrix  of  the  material  is  given  by 

1  \  ^  /  1 
COSV4  =1  cos(AV4)  —sin  (A  1/4) 

—i^esiniffe  cost/4  J  \  — /.^sin(AV4)  cos(AV4) 

A  similar  expression  for  TE  mode  can  be  obtained  by  replacing  the  3fe  with  the  negative  admittance  — Here  the 
and  ?Ve  are,  respectively,  the  generalized  impedance  and  admittance  because  they  include  incidence  angle,  i.e.,  Eq.  (HJi  is 
valid  for  both  normal  and  oblique  incidence.  They  are  given  by 


TM: 

TE: 


kl  =  klexliy-jkl, 
kl  =  kley^x-^kl 


(5) 


where  (Ev,  and  are,  respectively,  the  effective  permittivities  and  permeabilities  of  the  equivalent  layer. 

As  shown  in  Eq.  ©,  the  total  phase  of  a  A-layer  system  equals  A  times  the  phase  t/4  of  the  single  layer;  whereas 
the  impedance  or  admittance  ‘We  is  independent  of  the  number  of  layers.  These  properties  imply  that  the  material 
parameters  retrieved  from  the  A  layers  of  unit  cells  are  the  same  as  those  retrieved  from  one  layer  of  unit  cells.  In  other 
words,  the  bulk  metamaterial  parameters  can  be  predicted  from  a  single  symmetric  unit  cell.  Note  that  the  characteristic 
matrix  of  a  N-layer  asymmetric  unit  cells  does  not  have  the  nice  properties  as  those  in  Eq.  ®.  As  a  consequence,  the 
retrieved  material  parameters  will  be  dependent  on  the  number  of  unit  cells  along  the  propagation  direction.  Hence,  the 
bulk  material  parameters  cannot  be  predicted  from  one  unit-cell  layer.  In  other  words,  a  length-independent  description 
cannot  be  achieved  for  asymmetric  unit  cells. 


3.  Graphical  retrieval  method 

In  this  section,  we  will  introduce  a  graphical  retrieval  method  based  on  Herpin’s  theorem  and  the  previous  methodlfTTll. 
Erom  above  discussion  a  symmetric  unit  cell  can  be  represented  by  three  equivalent  layers  ABA.  Let  the  z-direction 
perpendicular  to  the  plane  of  the  layers.  Assume  each  layer  is  orthorhombic  anisotropic  with  the  principal  axes  parallel 
in  each  direction  and  the  effective  permittivity  and  permeability  tensors  are  given  by  Eq.  ©.  Let  the  x-z  plane  be  the 
plane  of  incidence.  The  characteristic  matrix  of  TM  mode  is  defined  as 


=  M  I 


Mil  Mi2\ 

M21  M22 J  J  ’ 


Ei 


(6) 


where  ‘i’  refers  to  the  input;  and  ‘o’  refers  to  the  output.  By  matching  boundary  condition,  the  relationship  between  the 
scattering  and  characteristic  matrices  of  a  unit  cell  is  given  by 


Mil  =  M22  = 


(1+52i)(1-5i2)+5ii522 


Mi2  = 
M21  = 


25ii 

(1+52i)(1+5i2)-5ii522 


(7) 


^  r 
Wx  - 


2Sn% 
{l-S2x){\-Sn)-SnS22 


where  2^i  and  2fo  are,  respectively,  the  generalized  input  and  output  impedance  whose  values  depend  on  the  incidence 
angle  and  background  permittivity  and  permeability  at  the  input  and  output  sides.  In  our  definition,  Sij{i  =  j)  is  related 
to  the  transmission  and  Sij{i  ^  j)  is  related  to  the  reflection.  That  means,  for  the  TE  polarization  and  5'2i  are, 
respectively,  the  transmission  and  reflection  coefficients  of  the  electric  field;  whereas  for  the  TM  polarization  the 


transmission  and  reflection  coefficients  of  the  electric  held  are  given  by 


£2^1: 


■^ii  and  — 5'2i,  respectively.  Note  that 


as  long  as  the  layer  has  inversion  symmetry.  Mu  =  M22  regardless  the  input  and  output  impedances.  For  scattering 
matrix,  5ii  =  S22  only  if  ^  In  addition  if  Mu  =  M22,  we  then  also  have  5i2  =  ^21.  From  Eq.  Q,  the  dispersion 

relation  and  the  effective  impedance  of  the  unit  cell  ABA  are  given  by 


cos{Ked)  =M\  \  =  cos{ki,di)cos{k2^d2)  —  T]^ sin{ki,di)sm{k2zd2) , 

■^21  +  77  +  cos(ki;^/i)sin(k2jC/2)  —  sin(k2zC/2) 


Mi2 


^2  sin(ki;t/i )  cos{k2zd2)  - 
sin(kij.c/i)cos(k2,c/2) - 


(8) 


77  +  cos(ki,^/i)sin(k2jC/2)  +  17  sm{k2^d2)  ’ 


where  the  phase  xj/e  =  Kgd  and  the  Ke  is  the  effective  propagation  constant.  The  subscript  1  and  2  refer  to  the  layer  A 
and  B,  respectively.  The  SI  is  the  generalized  impedance  of  the  equivalent  layer  A.  The  z-component  wave  vector  kij 
and  k2z  are  determined  from  the  dispersion  relations  in  Eq.  Q.  The  di  is  twice  the  thickness  of  layer  A;  and  the  d2  is  the 
thickness  of  layer  B.  The  period  of  the  symmetric  unit  cell  is  c/  =  c/i  +  d2-  Other  parameters  in  Eq.  ®  are  given  by 


±  _  1  /  ^2x^\z  ^  ^\x^2z 

2  \  ^\x  ^2z  ^2x  ^Iz 


B2x^\z  ^  BXx^2z 
BIx^2z  B2x^Iz 


for  TM , 
forTE. 


(9) 


Although  some  MMs  are  mesoscopic  media,  nevertheless  we  restrict  our  consideration  to  those  satisfying  d  after 
tedious  derivation,  Eq.  (H))  can  be  simplified  to 


El 

kl 

El 


-yr^x ' 


£xkl 

2^1  =  ^ 
£x 

kx 

kl£x£z 

Bxkl 

Bx 

kl 

'Pz^l  ’ 

klBxBz 

(10) 


Here,  the  ej  and  jj.j  {j  =  x,y,z)  are  the  bulk  values  of  the  effective  permittivities  and  permeabilities,  respectively.  They 
are  related  to  the  material  parameters  of  the  equivalent  layers  A  and  B  through 


Bp  = 


d-X^Xp  d2£2p 

d  ’ 
dxBxp  +  d2B2p 


£7  = 


d  £u£2i 


dx£2z  +  d2£xz 
d  Bxz  B2z 


(11) 


l^z  = 


dxB2z  +  d2Bxz 


where  p  =  x,y.  Equation  (fTTI)  is  a  simplification  of  the  most  general  linear  case  treated  by  Lakhtakialflhl.  In  practice, 
the  material  parameters  of  the  equivalent  layers  are  unknown.  One  can  only  measure  scattering  parameters.  In  most 
experiments  combining  Eqs.  (I?])  and  (fToli.  we  obtained  the  retrieval  formula: 


where 


TM: 

TE: 


Ym  =  e.vMy  -  , 

YE  =  eyT[,-^X, 


£x£z  ’ 


X  =  eblib&m^d,  Y,n  =  —Scos^9,  Ye  =  —Scos^9, 


Ym  =  Ye  = 


1 

kod 


£b 


2mn  ±  cos 


lib 


\-S: 


21  ■ 


25ii 


n  ^  2 


(1+521)^-52^ 
in  =  0,  ±1 ,  ±2,  •  •  • , 


(12) 


(13) 

(14) 


where  9  is  the  incidence  angle.  The  and  Pb  are,  respectively,  the  background  relative  permittivity  and  permeability. 
The  retrieval  formulas  in  Eq.  (fT^  provide  four  straight  lines  {Ym,  Y,„,  Ye,  Ye  vs.  X),  two  for  each  polarization.  The  Ym 
and  Ye  represent  the  dispersion  lines  for  TM  and  TE  polarizations,  respectively.  The  Y^  and  Ye  are  the  corresponding 
impedance  and  admittance  lines.  These  straight  lines  are  easy  to  implement  experimentally.  After  measuring  the 
scattering  parameters  at  several  incidence  angles  and  plotting  the  data  according  to  Eqs.  Call-lO,  use  linear  regression 
technique  to  calculate  the  slopes  and  Y-intercepts  of  the  four  lines.  Erom  the  slopes  and  Y-intercepts,  the  six  effective 
material  parameters,  Ej  and  Jij  {j  =  x,y,z),  can  be  retrieved  simultaneously.  Let  Y^  and  Yj^  represent  the  Y-intercepts  of 
the  two  lines  in  TM  polarization,  Sm  and  5m  the  corresponding  slopes  of  the  lines;  whereas  Y^,  Y^,  Se,  and  Se  are  the 
corresponding  quantities  for  TE  polarization.  Thus, 


TM: 

TE: 


± 


z  = 


± 


Y  = 


± 


liy 


lix 


Y  ’ 


(15) 


where  the  ±  sign  in  Eq.  (fTsTi  can  be  fixed  by  requiring  the  imaginary  part  of  refractive  index  greater  than  zero,  i.e. 
3(«m)  >  0,  3(ne)  >  0,  and  the  real  part  of  impedance  and  admittance  greater  than  zero,  i.e.  91  (Z)  >  0,  and  91(7)  >0  for 
passive  media ifTTl. 


4.  Resolving  phase  branch 

There  are  two  issues  regarding  the  branch  determination  in  Eq.  (fT4li.  the  ±  sign  in  front  of  the  inverse  cosine  and 
the  phase  branch  m.  The  inverse  cosine  cos^*(  )  represents  the  principal  value:  0  <  cos^*(  )  <  tt.  Knowing  only  the 
cosine  value,  the  phase  angle  cannot  be  determined.  One  need  the  information  of  sine,  which  can  be  obtained  from  the 
imaginary  part  of  the  cosine  as  following: 


A  = 


1  —  52 
i  ‘^21 


+  5?, 


25ii 


=  cos  ^  =  cos  ( ^,-  +  /^, )  =  cos  cosh  !  sin  sinh  (pi , 


(16) 


where  (pr  and  (p,  are  real.  Eor  passive  medium,  (p,-  >  0,  and  thus  sinh  0,  >  0.  Therefore,  the  ±  sign  in  front  of  the  inverse 
cosine  in  Eq.  (fT4li  can  be  resolved  from  the  imaginary  part  of  A: 


_f  cos-^A  if3(A)<0 

^  (2;r  — cos^'A  if3(A)>0' 


(17) 


After  determining  the  ±  sign,  next  step  is  to  resolve  the  correct  phase  branch  m.  This  part  can  be  very  confusing  when 
negative  refractive  index  might  be  involved.  Notice  that  among  the  four  retrieval  lines  in  Eq.  (fT2]).  only  the  dispersion 


lines,  Ym  and  Ye  ~  X,  are  functions  of  the  phase  branch  m,  the  impedance  and  admittance  lines  are  independent  of 
m.  Employing  this  property,  we  provide  three  methods  that  might  help  to  resolve  the  correct  phase  branch  m. 


Method- 1,  Using  e\  and  For  the  TM  polarization  in  Eq.  (fT^.  the  ej  can  be  obtained  either  from  the  intercepts 

Sm 

as  or  from  the  slopes  as  =  — ,  and  the  two  results  should  be  close.  When  the  branch  is  wrong,  the  two 

yO 

results  can  be  significantly  different.  Similarly  for  the  TE  polarization,  the  J[^  can  be  obtained  from  either  ^  or 


_ j  Se 

lit  =  — .  Hence,  the  correct  branch  m  is  the  one  that  minimizes  the  absolute  value  of  the  differences,  i.e. 

^5. 


TM: 

TE: 


mm 


mm 


S,n 

Y^{m)  SE{m) 


Y9 


m  =  0,  il ,  it2, 
OT  =  0,  il ,  it2,  • 


(18) 


Method-2,  Using  rr-  and  =- :  Alternatively,  we  can  use  ^  for  TM  and  for  TE  as  criteria  to  determine  the  correct 


phase  branch,  i.e. 


TM: 

TE: 


min  {\YM{m)Sm-Yl^,SM{in)\  :  m  =  0,  ±1,  ±2,  •  •  •  }  , 
min{\^Y^{m)Se —  Y^SE{m)\  :  m  =  0,±1, ±2, •  •  •  }  . 


(19) 


Method-3,  Using  EypLy  and  The  third  method  is  to  use  Ezpty  for  TM  and  EyH.  for  TE  as  criteria  to  select  the 

correct  phase  branch.  Then,  the  algorithm  becomes 


TM: 

TE: 


mm 


mm 


Suim) 

yHH 


yO 
^  m 

Sm 

vO 


Ssin 


:  OT  =  0,itl,±2, 
:  OT  =  0,itl,±2, 


(20) 


In  above  three  methods,  the  branch  number  m  predicted  by  the  first  two  methods  is  always  consistent  for  all  the 
frequencies  in  the  several  examples  we  tested.  The  third  method  predicts  the  same  result  as  the  first  two  for  most  of  the 
frequencies,  but  sometimes  it  can  be  different  by  ±  1  at  the  edge  of  phase  transition  frequencies  in  resonant  regimes.  Note 
that  before  applying  Eqs.  (fTSll.  (fTOb.  or  (l20l)  to  resolve  the  correct  phase  branch,  the  ±  sign  in  front  of  the  inverse  cosine  in 
Eq.  (fT4l)  must  be  determined  first.  From  the  several  examples  we  tested,  using  above  branch-resolving  techniques,  there 
is  no  need  to  recourse  adjacent  frequencies  in  determining  the  correct  phase  branch.  Recoursing  adjacent  frequencies 
can  become  confusing  when  both  positive  and  negative  refractive  indices  are  present  in  the  same  frequency  band. 


5.  Discussion 


In  this  section,  we  provide  an  example  to  show  how  to  implement  the  graphical  method  from  the  scattering  parameters. 
Typically  the  background  is  lossless  and  thus,  the  X  variable  in  Eq.  (fT2l)  is  real;  whereas  the  slops  and  Y-intercepts  are 
usually  complex.  The  real  and  imaginary  parts  of  the  retrieval  lines  should  be  plotted  separately  as  shown  in  Fig.|2]which 
was  calculated  from  the  scattering  matrix.  Drude  model  is  used  for  the  effective  material  parameters  of  the  equivalent 
layers  A  and  B, 


e  =  1  -  ■ 


r 

J  ep 


r 

Jmp 


^  P-flr+m' 

where  fgp  =  30THz,  /„,p  =  20THz,  and  y  =  3THz.  The  Ex  and  jXx  are  described  by  Eq.  d2Tli  with  resonances  at 
fer  =  20THz  and  f„tr  =  25THz  for  the  layer  A  and  at  fer  =  35THz  and  f„tr  =  37THz  for  the  layer  B.  The  other 


(a)  TM  polarization:  dispersion  lines  (left  panels)  and  impedance  (b)  TE  polarization:  dispersion  lines  (left  panels)  and  admittance 

lines  (right  panels).  lines  (right  panels). 


Fig.  2.  Retrieval  lines  at  the  frequency  30THz  calculated  from  the  scattering  matrix  of  a 
single  layer  of  unit  cells  at  seven  incidence  angles  (denoted  as  circles)  uniformly  distributed 
from  0  to  30  degree.  The  horizontal  axis  is  sin^S.  Background  medium  is  vacuum.  Top 
panels:  real  part.  Bottom  panels:  imaginary  part. 


parameters  for  A:  Ey  =  Ex  —  0.3,  e,  =  +  2,  fly  =  fix  —  0.5,  and  =  1.  The  other  parameters  for  B:  Ey  =  Ex  —  0.8, 

Ey  =  Ex  —  0.5,  jiy  =  IJ.X  +  0.2,  and  jj,,  =  fix  —  0.6.  The  thickness  is  240nm  for  the  layer  A  and  320nm  for  the  layer  B. 
Thus,  the  period  of  the  unit  cell  (ABA)  is  800  nm.  Figure^shows  the  correct  phase  branch  m  predicted  for  each  frequency 
in  the  regime  of  interest  for  the  thickness  of  one  (a)  and  six  (b)  periods.  In  our  example  since  d  <^X,  for  one  unit-cell 
thickness,  most  of  the  frequencies  are  within  the  fundamental  branch  m  =  0  except  for  the  regime  of  negative  index  of 
refraction  [see  Fig.Ua)]  where  m  =  —1.  For  the  six-period  thickness,  the  phase  branch  jumps  to  m  =  1  in  the  frequency 
range  18  ^  20THz  due  to  the  high  valves  of  the  positive  refractive  index  in  this  regime[see  Fig.Ua)]. 

Shown  in  Fig.  |4]  are  the  effective  index  of  refraction,  impedance,  and  admittance  retrieved  from  one  unit  cell.  These 
values  are  the  same  as  those  retrieved  from  the  six  unit  cells.  In  this  example,  both  TM  and  TE  polarizations  contain  two 
frequency  bands  where  the  effective  index  of  refraction  is  negative.  As  illustrated  in  Fig.  [3  these  two  negative  bands 
are  accurately  captured  by  the  branch-resolving  techniques  provided  in  Sec.  01  The  six  retrieved  effective  permittivities 
and  permeabilities  are  shown  in  Figs.|5]and|6]in  the  blue-solid  curves,  which  are  obtained  from  the  graphical -retrieval 
and  phase-unwrapping  scheme,  i.e.  Eqs.  (ITSll.  (ITtI).  and  (l20l).  Eor  comparison,  we  also  show  the  effective  permittivities 
and  permeabilities  calculated  from  Eq.  (ITTli  in  Eigs.|5]and|6]in  the  green-circle  curves.  The  retrieved  effective  material 
parameters  agree  very  well  with  the  results  calculated  from  the  effective  medium  theory.  These  values  represent  the  bulk 
values  of  the  material  parameters  since  they  are  extracted  from  the  symmetric  unit  cells.  We  have  successfully  applied 
the  graphical  method  and  phase-unwrapping  techniques  to  our  recent  experimentlflTl.  Erom  the  experimental  point  of 
view,  the  straight-line  graphical  method  is  more  accurate  than  the  methods  using  single  data  point.  The  linear  regression 
technique  is  based  on  collective  data  points  measured  at  several  incidence  angles.  It  has  an  averaging  effect  and  thus 
reduces  the  uncertainty  of  measurements.  Using  several  incidence  angles  can  also  help  to  resolve  the  phase  ambiguity. 
After  determining  the  effective  parameters  of  the  bulk  material,  Eq.  (ITTI)  can  be  used  to  recovery  the  effective  parameters 
of  the  equivalent  layers  if  we  are  interested  in.  This  extra  benefit  may  help  to  improve  the  unit-cell  design. 

When  the  tensors  in  Eq.  O  rotate  about  the  z-axis  without  rotating  the  coordinate  system,  the  off-diagonal  elements 
in  the  x-y  plane  will  not  be  zero  {Exy  ^  0  and  fixy  ^  0),  and  thus  the  x  and  y  components  of  the  electromagnetic  fields 
will  be  coupled.  The  dispersion  relations  in  Eq.  (|5]l  will  no  longer  be  valid.  Without  the  proper  modification,  the  current 
parameter  retrieval  scheme  cannot  be  applied  to  this  scenario.  The  original  classification  of  the  left-handed  and  right- 


TM:  Number  of  Period  =  1  TM:  Number  of  Period  =  6 


(a)  One-period  thickness.  (b)  Six-period  thickness. 

Fig.  3.  Branch  number  predicted  by  the  method-3  [Eq.  (doll]  vs.  frequency.  Since  d  X, 
most  of  the  frequencies  are  in  the  fundamental  branch  m  =  0.  The  m  =  —  I  indicates  the 
frequencies  of  the  negative  refractive  index  [n  <  0,  see  Fig.|4(a)|.  The  m  =  1  corresponds 
to  the  frequencies  of  the  high  valves  of  the  positive  refractive  index  [see  Fig.  |4(a)|.  Top: 

TM  polarization.  Bottom;  TE  polarization. 


(a)  Retrieved  refractive  index  for  TM  polarization  (left  panels)  and  (b)  Left  panels:  Retrieved  impedance  for  TM  polarization.  Right 

for  TE  polarization  (right  panels).  panels:  Retrieved  admittance  for  TE  polarization. 


Fig.  4.  Effective  refractive  index,  impedance,  and  admittance  retrieved  from  a  single  unit¬ 
cell  layer.  Top  panels:  real  part.  Bottom  panels:  imaginary  part. 


handed  electromagnetic  materials  can  cause  confusion  with  chiral  materials  which  are  important  class  of  electromagnetic 
materials ifTsl.  To  avoid  such  confusion,  we  can  use  phase-velocity  characteristics.  When  the  phase  velocity  and  the 
time-averaged  Poynting  vector  form  an  acute  angle,  the  medium  is  positive-phase-velocity  characteristics.  When  the 
phase  velocity  and  the  time-averaged  Poynting  vector  form  an  obtuse  angle,  the  medium  is  negative-phase-velocity 
characteristics. 


(a)  The  effective  (left  panels)  and  (right  panels).  (b)  The  effective  Ey  (left  panels)  and  fly  (right  panels). 

Fig.  5.  In-plane  values  of  the  retrieved  material  parameters.  Upper  panels:  real  parts.  Lower 
panels:  imaginary  parts. 


Fig.  6.  The  effective  (right  panels)  and  (left  panels).  Upper  panels:  real  parts.  Lower 
panels:  imaginary  parts. 


6.  Conclusions 

In  conclusions,  we  have  shown  that  symmetric  unit  cells  can  often  be  constructed  in  metamaterials  regardless  the  number 
of  sub-layers  and  the  complexity  of  each  sub-layers  in  the  original  unit  cells.  The  graphical  retrieval  method  and  phase 
unwrapping  techniques  presented  here  might  be  a  useful  tool  for  metamaterial  designs. 
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